Analytical and Numerical Modeling of Solute Intrusion,Recovery, and Rebound in Fractured Bedrock     

Ranjeet M. Nagare  Young-Jin Park  Trevor Butterfield  Chad Belenky  Sheldon Scyrup 《Ground water》2020,58(1):56-69

Contaminated groundwater in fractured bedrock can expose ecosystems to undesired levels of risk for extended periods due to prolonged back-diffusion from rock matrix to permeable fractures. Therefore, it is key to characterize the diffusive mass loading (intrusion) of contaminants into the rock matrix for successful management of contaminated bedrock sites. Even the most detailed site characterization techniques often fail to delineate contamination in rock matrix. This study presents a set of analytical solutions to estimate diffusive mass intrusion into matrix blocks, it is recovered by pumping and concentration rebound when pumping ceases. The analytical models were validated by comparing the results with (1) numerical model results using the same model parameters and (2) observed chloride mass recovery, rebound concentration, and concentration in pumped groundwater at a highly fractured bedrock site in Alberta, Canada. It is also demonstrated that the analytical solutions can be used to estimate the total mass stored in the fractured bedrock prior to any remediation thereby providing insights into site contamination history. The predictive results of the analytical models clearly show that successful remediation by pumping depends largely on diffusive intrusion period. The results of initial mass from the analytical model was used to successfully calibrate a three-dimensional discrete fracture network numerical model further highlighting the utility of the simple analytical solutions in supplementing the more detailed site numerical modeling. Overall, the study shows the utility of simple analytical methods to support long-term management of a contaminated fractured bedrock site including site investigations and complex numerical modeling.  相似文献   

Gondwana-Derived Terranes Structural Mapping Using PALSAR Remote Sensing Data   总被引：1，自引：0，他引：1  

Amin Beiranvand Pour  Mazlan Hashim  Yongcheol Park 《Journal of the Indian Society of Remote Sensing》2018,46(2):249-262

Gondwana-derived terranes are now separated by major faults or suture zone, which are rarely simple and easily recognizable lineaments. Different association of ore mineral systems such as deposits of sediment-hosted/orogenic gold and granite-related minerals is discovered in collision and subduction zones of the Gondwana-derived terranes. They are associated with large-scale, terrane-bounding fault systems and broad areas of deformation. Mineralization mostly associated with structurally-controlled complex lodes, veins, sheeted veins and veinlets in diverse orientations. Recognizing the structural significance of lineaments and curvilinear is very difficult in tropical, arid and Antarctic regions due to environmental obstacles. Remote sensing data could be used to detect geological structures associated with suture zones between Gondwana-derived terranes especially for large inaccessible regions where fieldwork is limited or nonexistent. In this investigation, the Phased Array type L-band Synthetic Aperture Radar (PALSAR) satellite remote sensing data were used to map major geological structures in the Bentong-Raub Suture Zone of Peninsular Malaysia. It is one of the major structural zones in Sundaland, Southeast Asia, which forms the boundary between the Gondwana-derived the Sibumasu terrane and Sukhothai-Indochina arc. Structural features associated with sediment-hosted/orogenic gold deposits in the Central Gold Belt and tin mineralization in S-type granites in the Main Range were investigated using PALSAR data and comprehensive fieldwork. Results indicate that main faults strike along 340° to 350° that are intersected by many shear or lateral fault zones are high potential zone for gold mineralization in the Central Gold Belt. Hydrothermally alteration mineral zones and cataclastic rocks are also the other indicators of gold mineralized veins in the gold belt. High potential zones for tin mineralization are N–S strike-slip faults, fault zones and shear zones trending E–W, NE–SW and WNW–ESE in dissected crystalline granitic rocks that are associated with hydrothermal alteration zones in the Western Tin Belt. Lineament analysis using PALSAR satellite remote sensing data is a useful tool for mapping major geological structural features and detection of the boundary between the Gondwana-derived terranes and detailed structural analysis of fault systems and deformation with high potential for a variety of mineral resources, especially in tropical, arid and Antarctic regions.  相似文献   

Day-Night and Gastrointestinal Changes in Winter Feeding Intensity and Prey Composition of Apogon lineatus Inhabiting the Southeastern Waters of Korea     

Park  Joo Myun  Jung  Hae Kun  Lee  Chung Il  Huh  Sung-Hoi 《Ocean Science Journal》2019,54(1):141-146

Ocean Science Journal - Diel and gastrointestinal changes in relation to the feeding intensity and prey composition of Apogon lineatus were investigated through the analysis of the diet contents of...  相似文献   

Comparison of abstraction scenarios simulated by SWAT and SWAT-MODFLOW     

Eugenio Molina-Navarro  Ryan T. Bailey  Hans Estrup Andersen  Hans Thodsen  Anders Nielsen  Seonggyu Park 《水文科学杂志》2019,64(4):434-454

In hydrological modelling of catchments, wherein streams are groundwater-fed, an accurate representation of groundwater processes and their interaction with surface water is crucial. With this purpose, a coupled model was recently developed linking SWAT (Soil and Water Assessment Tool) with the fully-distributed groundwater model MODFLOW (Modular Groundwater Flow). In this study, SWAT and SWAT-MODFLOW were applied to a Danish groundwater-dominant catchment, simulating groundwater abstraction scenarios and assessing the benefits and drawbacks of SWAT-MODFLOW. Both models demonstrated good performance. However, SWAT-MODFLOW provided more realistic outputs when simulating abstraction: the decrease in streamflow was similar to the volume of water abstracted, while in SWAT the impact was negligible. SWAT also showed impacts on streamflow only when abstractions were taken from the shallow aquifer, not from the deep aquifer. Overall, SWAT-MODFLOW demonstrated wider possibilities for groundwater analysis, providing more insights than SWAT in supporting decision making in relation to environmental assessment.  相似文献   

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock          下载免费PDF全文

J.‐W. Park  C. Park  J.‐W. Song  E.‐S. Park  J.‐J. Song 《国际地质力学数值与分析法杂志》2017,41(6):880-898

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

SHRIMP U–Pb ages of detrital zircons from the Early Cretaceous Nakdong Formation,South East Korea: Timing of initiation of the Gyeongsang Basin and its provenance          下载免费PDF全文

Tae‐Ho Lee  Kye‐Hun Park  Keewook Yi 《Island Arc》2018,27(5)

To constrain the depositional age of the lowermost Nakdong Formation in the Early Cretaceous Gyeongsang Basin, SHRIMP U–Pb age determination was carried out on zircon separates. The U–Pb compositions of detrital zircons from the Nakdong Formation yield a wide range of ages from the Archean to the Cretaceous but show a marked contrast in age distribution according to the geographical locations within the basin. The provenance of the southern Nakdong Formation is dominantly the surrounding Yeongnam Massif, which is composed of Paleoproterozoic metamorphic rocks and Triassic to Jurassic plutonic rocks, whereas the central to northern Nakdong Formation records significant sediment derivation from the Okcheon Metamorphic Belt, which is distributed to the northwest, in addition to the contribution from the Yeongnam Massif. It is suggested that the maximum depositional age of the Nakdong Formation is ca 127 Ma, based on its youngest detrital zircon age population. The onset of its deposition at 127 Ma coincided with the tectonic inversion in East Asia from a compressional to an extensional geodynamic setting, probably due to the contemporaneous change in the drift direction of the Izanagi Plate and its subsequent oblique subduction.  相似文献   

Effects of fluid displacement pattern on complex electrical impedance in Berea sandstone over frequency range 104–106 Hz          下载免费PDF全文

Yi Zhang  Hyuck Park  Osamu Nishizawa  Tamotsu Kiyama  Yu Liu  Kwangseok Chae  Ziqiu Xue 《Geophysical Prospecting》2017,65(4):1053-1070

To better understand the effect of fluid distribution on the electric response of rocks saturated with oil and brine, we conducted experimental studies on the complex electrical impedance in a Berea sandstone, together with in situ acquisitions of oil distribution images employing a high‐resolution medical X‐ray computed tomography. We performed two tests of brine displacement by oil under high (10 MPa) and low (5 MPa) pressures, which were accompanied by fingering and stable displacement patterns, respectively. The measured complex impedance data were fitted to the Cole model to obtain the resistance, capacitance, peak frequency of the imaginary impedance, and the exponent α of the rock–fluid system. With increasing oil saturation, the resistance showed an increasing trend, whereas the other three parameters decreased. The fingering displacement exhibited lower resistance and capacitance than the stable displacement. The analysis of the resistance changes using a simple parallel connection model indicates that there are more components of residual brine in parallel connections in the fingering pattern than in the stable displacement pattern at the same saturation. We also interpreted the normalised changes in the capacitance (or apparent dielectric constant) with respect to the oil saturation via an analysis of the shape factor of fluid distribution based on the Maxwell–Wagner–Brugermann–Hanai model. The changes in the shape factor suggest that the pinch‐off of the brine in parallel connection by the oil is a dominant mechanism reducing the capacitance. In the stable displacement, most of the connections in the brine phase are immediately pinched off by oil displacement front at a local oil saturation of 65%. Conversely, in the fingering displacement, there is a transition from the bulk or layered brine to the pinched‐off at a local oil saturation below 60%. The analyses indicate that the difference in the fluid distribution under different fluid conditions is responsible for the non‐Archie behaviour.  相似文献   

Relation of CME Speed and Magnetic Helicity in CME Source Regions on the Sun during the Early Phase of Solar Cycles 23 and 24     

R.-S. Kim  S.-H. Park  S. Jang  K.-S. Cho  B. S. Lee 《Solar physics》2017,292(4):66

To investigate the relations between coronal mass ejection (CME) speed and magnetic field properties measured in the photospheric surface of CME source regions, we selected 22 disk CMEs in the rising and early maximum phases of the current Solar Cycle 24. For the CME speed, we used two-dimensional (2D) projected speed observed by the Large Angle and Spectroscopic Coronagraph onboard the Solar and Heliospheric Observatory (SOHO/LASCO), as well as a 3D speed calculated from the triangulation method using multi-point observations. Two magnetic parameters of CME source regions were considered: the average of magnetic helicity injection rate and the total unsigned magnetic flux. We then classified the selected CMEs into two groups, showing: i) a monotonically increasing pattern with one sign of helicity (group A: 16 CMEs) and ii) a pattern of significant helicity injection followed by its sign reversal (group B: 6 CMEs). We found that: 1) 3D speed generally shows better correlations with the magnetic parameters than the 2D speed for 22 CME events in Solar Cycle 24; 2) 2D speed and the magnetic parameters of 22 CME events in this solar cycle have lower values than those of 47 CME events in Solar Cycle 23; 3) all events of group B in Solar Cycle 24 occur only after the beginning of the maximum phase, a trend well consistent with that shown in Solar Cycle 23; 4) the 2D speed and the helicity parameter of group B events continue to increase in the declining phase of Solar Cycle 23, while those of group A events abruptly decrease in the same period. Our results indicate that the two CME groups have a different tendency in the solar cycle variations of CME speed and the helicity parameters. Active regions that show a complex helicity evolution pattern tend to appear in the maximum and declining phases, while active regions with a relatively simple helicity evolution pattern appear throughout the whole solar cycle.  相似文献   

Beta-�� distribution and its application to hydrologic events     

Md. Sharwar Murshed  Sangil Kim  Jeong-Soo Park 《Stochastic Environmental Research and Risk Assessment (SERRA)》2011,25(7):897-911

The beta-κ distribution is a distinct case of the generalized beta distribution of the second kind. In previous studies, beta-p and beta-κ distributions have played important roles in representing extreme events, and thus, the present paper uses the beta-κ distribution. Further, this paper uses the method of moments and the method of L-moments to estimate the parameters from the beta-κ distribution, and to demonstrate the performance of the proposed model, the paper presents a simulation study using three estimation methods (including the maximum likelihood estimation method) and beta-κ and non beta-κ samples. In addition, this paper evaluates the performance of the beta-κ distribution by employing two widely used extreme value distributions (i.e., the GEV and Gumbel distributions) and two sets of actual data on extreme events.  相似文献   

Modeling a major source of economic resilience to disasters: recapturing lost production     

JiYoung Park  Joongkoo Cho  Adam Rose 《Natural Hazards》2011,58(1):163-182

Terrorist attacks and natural disasters have potentially severe economic consequences in terms of property damage and business interruption. However, experience from the September 11 World Trade Center attack and other disasters indicates that the economy has a great deal of resilience. This refers to the ability to dampen the maximum potential economic output (business interruption) loss. One of the most prominent sources of resilience is the ability of businesses to reschedule, or recapture, lost production after the event. Although there have been applications of a fixed parameter recapture factor for each of several aggregated sectors of the economy, there has been little formal analysis of this resilience action. This study offers a theoretic framework for analyzing production rescheduling. It distinguishes the major conditions influencing two aspects that have previously been neglected: (1) the maximum time span over which the rescheduling can take place and (2) the likely decline of the maximum recapture as the business interruption increases. We divide the relevant time path into two periods after recovery. One is a function of a recaptured output path after recovery to the status of normal production. The other is a function for the maximum recaptured production, based on the recaptured output path. The recaptured output path function is assumed to follow a normal distribution function, and hence, total recaptured output follows the cumulative normal distribution function over time after productive capacity is restored. Also, we develop a new cumulative normal distribution function for interruption time duration, which is symmetric with respect to the output axis. This recapture function has unknown parameters. Empirical data on the recaptured amounts following an actual disaster can be used to estimate the parameters of this function using simulation methods.  相似文献